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Heavy flavours represent a challenge for lattice QCD. We discuss it in very general terms. 
We give an idea of the significant recent progress which opens up good perspectives for high 
precision first principles QCD computations for fiavour physics. 

1 Overview 

In this talk I want to convey an idea about the perspectives for precise lattice QCD computations. 
Some emphasis is put on the field of flavour physics, where lattice QCD seems to be needed 
the most in the quantitative interpretation of present and future experiments (in particular this 
is the reason for Sect. 14. 3p . I do not want to discuss the latest numbers but will focus on the 
principle, the challenges and the perspectives. I mainly address those in the audience who know 
little about the field and therefore take a bit of a bird's-eye view of the field. So I do not 
hesitate to also list trivial facts (e.g. Sect. 14. 2p and my personal opinion on what is relevant for 
the future. 

2 The principle 

We want to find answers for properties of QCD in the non-perturbative domain, where the usual 
power series expansion in the coupling constant fails completely. This is done through a 

Monte Carlo (MC) evaluation 

of the Euclidean path integral 

after a discretization of space-time on a lattice with spacing a in all 3-1-1 dimensions. 

The last of these items seems like a bad mutilation of the theory, but we emphasize that dis- 
cretization is just one form of regularization, i.e. the taming of the infamous UV divergences 
of continuum path integrals. In fact discretization is the only known, complete regularization 
of the theory: complete in the sense that it defines the theory mathematically also beyond an 
expansion in the coupling. 



In lattice QCD the Euclidean path integral is "solved" numerically. Many observables can 
be computed directly in Euclidean space-time, but even more cannot be computed (unless the 
complete dependence on time is known). Examples are scattering cross sections. 

What at first sight looks like the (classical) coupling constant, go, and quark masses, mo,i, 
have no direct physical meaning in a quantum field theory such as QCD. They are bare (actu- 
ally divergent) quantities which need renormalization. In the general (non-perturbative) case, 
renormalization just means that the bare parameters are replaced by observable quantities. It is 
convenient to parametrize the theory as follows. The MC evaluation of the path integral yields 
dimensionless hadron masses (remember that a is the lattice spacing), 

am)i-'i = Mr(50,{/xr^'^}) (1) 

as a function of the bare, dimensionless, parameters in the Lagrangian of the theory, go, jj'^^^'^^ = 
^"^oT'^'^ , i = u,d, . . .. We consider the theory with Nf quarks and require A^f ratios of hadron 
masses to agree with Nature. This may be read as fixing the bare dimensionless quark masses 
/Xj to the correct values at a given go. Taking now in addition the dimensionfull masse(s) ruh 
from experiment, yields the lattice spacing 

had 

a = ^ = a{go) , (2) 
rrih 

for any value of the bare coupling go- At this point the theory is parametrized in terms of 
physical observables namely A'^f + 1 hadron masses. A new observable, e.g. the B-meson decay 
constant, F^, relevant for flavour physics, is now a function 

i^B = i^B(o,mp,m^, . . .) . (3) 

With the typical approximation of neglecting the u-d quarkmass difference as well as the influ- 
ence of the top-quark on low-energy QCD, a useful set to be chosen in eq. ([3|) is 



= 938.3 MeV , = 139.6 MeV , rriK = 493.7 MeV , mo = 1896 MeV , me = 5279 MeV . 



We observe that this experimental input is very precise. Using it as the starting point, the 
running renormalized quark masses and coupling can be computed (at high momentum scale) 
as well as other quantities; in particular matrix elements needed for flavour physics, such as F^. 

Of course the continuum limit a ^ has to be taken to arrive at the physical quantity, 
-Fb(0, ?7ip, rnvr, . . .). The only assumption is that this limit exists. (In practise, it is taken by a 
numerical extrapolation of a few data points at different a.) There is a large amount of evidence, 
perturbative and non-perturbative and also from the investigation of critical phenomena, that 
given a local action (that means a discretization with exponentially localized interactions) the 
continuum limit does exist. The perturbative renormalization group self-consistently predicts 
that the continuum limit is at go = 0: a(0) = 0. 

2.1 Some results 

In recent years more and more results have appeared with rather small quote d er ror bars. We 
refer the reader to reviews by M. Delia Morte and A. Jiittner at "Lattice 2007"[Il2land, focusing 
more on recent results, the one of E. Gamiz at "Lattice 2008" l^. In Tabled] we collected a few 
quantities, interesting for flavour physics, for which errors around 1 % are quote di 



"With <^B = ciFb, the equivalent dimension-less form is 



— — -Rb(Mp: — ) • 



and -Fb ~ mpRsiO, j^, ■ ■ ■) is the continuum prediction. 



m^{3 GeV) = 0.986(10) GeV HPQCdI^I 

mp(mb) = 4.20(4) GeV HPQCdEI 

^ = ^§^^f- = 1.211(38)(24) FNAL/MILCEI 

Fb, = 243(11) MeV FNAL/MILcEl 

Fd, = 241(3) MeV HPQCdEI 



Table 1: Results discussed in the review of E. Gamiz' 



m 



Taken at face value, the last line in the table is in disagreement with what the CLEO 
experiment finds; see the talk of Philip Rubin at this conference. Should we conclude in favour 
of new physics^^, or is this precocious? As a preparation for an answer to this question, we take 
a rough look at the machinery used in practise to obtain the results. 

2.2 The machinery 

At this point in time, the numbers with the smallest quoted error bars come from so-called 
rooted staggered fermion simulations (all those in Table [T] do). There are serious concerns 
about the validity of this formulation, see for example^. Since a local action equivalent to 
rooted staggered fermions is not known, the argume nts iii stifving the rooting trick represent an 
involved discussion based on additional assumptions. ^^^^^ 

Also NRQCD is used for the b-quark in the analysis which yields Table [TJ As an effective 
field theory it contains power divergences; this means terms that diverge as l/a"". In lattice 
NRQCD they are subtracted perturbatively. Uncancelled divergences of the form 

I :^oo (4) 

a ruf, a [log(a)]'^mb 

are then left behind when one works until /c'th order in perturbation theory. The continuum 
limit therefore does not exist and the estimate of the precision of extracted results is delicate. 
Other sophistications commonly employed include smearing, many parameter fits (in the case 
of staggered chiral perturbation theory) and Baysian fits. 

All in all the large amount of tricks employed in the computations compromises their direct- 
ness and "first principles status" . One could also say systematic errors might be underestimated. 
It is often emphasized that the methods are tested by a comparison to experimental numbers. 
In our opinion, this is appropriate for gaining a rough idea on the validity range of model calcu- 
lations, but not the way to establish a controlled theory computation. Before concluding that 
new physics is present, it appears very desirable to perform independent computations with an 
independent technology. In particular we insist on 

• a manifestly local discretization and 

• non-perturbative subtraction of power law divergences. 

Such computations are in progress. Before turning to them, let us understand why they were not 
carried out already a number of years ago. The basic reason is that lattice QCD is a considerable 
challenge. 

3 The challenge 

This becomes apparent by considering just the typical hadronic input eq. ([3|). There is a large 
range of scales between m-,^ ~ 140 MeV over mo = 2 GeV to me = 5 GeV. In addition, the 
ultraviolet cutoff of Auv = of the discretized theory has to be large compared to all physical 
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Figure 1: Performance of DD-HMC as a function of the inverse quark mass' — '. "Accelerated DD-HMC" refers 
(amongst others) to the use of deflation. Note that the cpu-time t for original HMC algorithm scales with 



n > 3. 



energy scales if the discretised theory is to be an approximation to a continuum. Also the linear 
extent of space time has to be restricted to a finite value L in a numerical treatment: there is 
an infrared cutoff L~^. Together the following constraints have to be satisfied. 

Air = L^^ < , . . . , rriD , me <C a^^ = Auv 
effects 0(e-^™-) 0{{am^^'^f) 

I I 
L>4/m^^6im a < l/(2mD) ~ 0.05 fm 

yielding L/a> 120. 

After the first line, we have discarded the scale of mesons or baryons with b-quarks. Thus 
these are not included in the already rather intimidating estimate of L/a > 120. Looking at 
these numbers, it appears unavoidable to separate the b-quark mass scale from the others before 
simulating the theory. We shall return to this in Sect. 14.31 

4 Perspectives 

Despite the considerable challenge that we are facing, there are good perspectives for meeting it 
soon. First, without being able to enter into the merited in-depth discussion, I mention that the 
condition L > 6 fm may be relaxed by choosing unphysically large pion masses combined with a 
subsequent extrapolation. A factor of 1.5 to 2 maybe gained this way. Such extr ap olat i ons are 
theoretically guided by chiral perturbation theorjA^, also including lattice spacing effectsl^^^^. 

4-1 Algorithms 

Second, it is not obvious what is the CPU-effort for a QCD simulation. Indeed, since the 
development of the first exact algorithm called Hybrid Monte Carlo (HMC)^!^, it has been 
improved consi derably | by multiple time-scale integration^I^IIS], the Hasenbusch trick of mass- 
preconditioning f^SEDI ii^Q of the domain decomposition method in QCD ^^^ 33 qj^^ 



''Even though we list only buzzwords here, it is impossible to do justice to all relevant developments and the 
interested reader is advised to consult reviews^^. 
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Figure 2: HQET strategy. The parameters in the HQET Lagrangian and fields are denoted by oj^. In the matching 
step they are directly determined from requiring ^^'^^^ {a, Li,mh) Li,mh) (other masses are 



supressed). In a second step they are used to compute $j ("^ = 0,1/2 = 2Li,mb) by an extrapolation in ai. 
Returning to finite a, allows to extract ujk{a2) at larger lattice spacings a2. These are sufficient for infinite volume 
computations. The r.h.s. shows the same procedure understood as a relation between renormalized quantities. 



deflation method combined with chronological inverters In parallel, parameters of the 
algorithms may be tuned^^Sl and it has been understood that instabilities due to the breaking 
of chiral symmetry become irrelevant in sufficiently large volumes'^. Also algorithms for the 
contri butions of single flavours (not mass-degenerate within a doublet) have been developed 
l ^^l^^ l It is further su ggested that fluctuations at small quark masses can be tamed by clever 



' ' '. it IS lurtner su ggested 
reweighting techniquesl^SEQ]^ 



Amongst all these developments, we emphasize the recent taming of the critical slowing 
down with the quark mass. It was achieved by M. Liischer through a combination of deflation 
with others of the tricks listed above. The impressive weak dependence of the execution time 
on the quark mass is illustrated in Fig. [1] taken from^^^l^ 



4-2 Machines 

It is well known that the speed of computers has been increasing continuously - at an exponential 
rate. Yet, it is illustrative to look at an (incomplete and personal) list of the machines available 
over the years. Tflop numbers have to be taken with care, but the overall message on the fast 
growth of available computer power is clear. Because of recent investments into high performance 
computing, the growth has been stronger than Moore's law. 



year 


machine 


speed/Tflops 


share for a typical 
LQCD collaboration 


1984 


Cyber205 


0.0001 


/lOO 


1994 


APEIOO 


0.0500 


/4 


2001 


APEIOOO 


0.5000 


/4 


2005 


apeNEXT 


2.0000 


/3 


2009 


PC-Cluster "Wilson" 


15.0000 


/I 


2009 


BG/P 


1000.0000 


/20(?) 



Table 2: An incomplete list (in particular between 1984 and 1994 there were a number of other machines) of 
computers on which we have done QCD computations. 



4-3 Strategies for the b-quark 



B-meson decays and mixing are an essential part of flavour physics. We liave seen above that 
their direct, "relativistic" , treatment is impossible, since the demand of L/a > 120 would be 
increased by yet another factor of 2-3. Conversely we may, however, also use the large mass 
of the b-quark to our advantage, working in an expansion in AqcD/mt,. There are kinematical 
limitations to this expansion: momenta of other particles also have to be small compared to the 
b-quark mass. This means for example that the larger part of phase space in a decay B ttIv 
is not accessible to HQET, the effective field theory which implements the expansion. In the 
analysis of the experimental data an apropriate cut has to be implemented. 

The most significant difficulty in a non-perturbative use of HQET is, however, of a different 
origin. We already mentioned power law divergences in Sect. 12.21 These appear in the effective 
theory and have to be subtracted non-perturbatively, which means the parameters of the effective 
theory have to be determined non-perturbatively, thus avoiding terms of the form eq. (jH). In 
contrast to QCD, new parameters appear at each order in the HQET expansion. Still, their 
determination follows the logic as described around eq. ([T|). There is only one twist to the general 
principle: instead of using a relatively large number of experimental inputs, which results in a 
loss of predictive power |, one may compute quantities in lattice QCD with a relativistic b-quark 
and use those as in put. T his step is denoted by non-perturbative matching. It was realized 
already a while ago^^^^^S] ^j^at such a strategy can be carried out by considering quantities 
defined in a finite volume of linear dimension around L = 0.2 . . . 0.5 fm = 0(1/Aqcd) in this 
initial step. Some technical obstacles had to be removed^^^IMIMIMI^ i^^^ j^ow the path appears 
paved for a precise implementation of the idea. The strategy is depicted in Fig. [2j 

An impression of the achievable precision may be gained from the pioneering computation 
fSSlj where for the first time l/?n.b corrections were taken into account. We warn the reader that 
this was a quenched computation, where sea-quark effects are ignored. It should be used to judge 
on the achievable precision, not for input into phenomenology. The computation focussed on 
the relation between the (spin-averaged) B-meson mass and the renormalization group invariant 
b-quark mass, -Mb. Table [3] displays the results for M^, at the lowest order in l/m-b. Units are 
the potential scale^^Zl ^.^ ~ 0.5 fm. The second column of the table lists the results in the static 
approximation, where terms of order A in Mb are included, while corrections of order A^/m^ 
relative to the leading term are omitted. Since the computation is carried out with different 
matching conditions, specified by an angle appearing in the boundary conditions for the 
quarks, their spread of about 3% can be taken as an indication for the size of the corrections. 
At the next order in l/rrib (columns 3-5), even more different matching conditions are explored 
(angles 61,02). The spread of results is now reduced to a level of below 1 part in 200, much 
smaller than the total error of the result @. We conclude that precision computations are possible 
in HQET. 

As a preparation for a discussion of a second method for implementing the scale separation, 
let us turn to the right part of Fig. [2l There the procedure is interpreted as a relation 

^.QCD(ao,Li,mb) = cI>HQi^T(ai,Li,mb) - ^"^'^^(ai, L^, r^b) = ^''^''^ {a^, L2,m^) 

large volume physics 

The fact that only observables appear in this chain, with the parameters in the Lagrangian of 
the theory eliminated, makes clear that the continuum limits ^ can be taken in each step. 
In practise, oq down to 0.01 fm is used while 02 ranges between 0.05 fm and 0.1 fm. 



"^For example it is hard to imagine how one could avoid to have the decay constant among the input instead 
of among the predictions. 

'^Such a statement can be made since errors for different entries in the table are signficantly correlated. 
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0.6% ^ total error 
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Table 3: 


Renormalization 


group invariant b-quark mass 


. A = Aqcd. 



The second method starts from the same idea and the same (or similar) quantities in L = Li. 
But then it remains in relativistic QCD, computing 

oALi) = lim lim l.^^^^ ^' — — ^ , (5) 

for each desired observable ^S*""^. Similarly one connects to the large volume in the last step 
I '^^ l '^^ l The advantage is that there are no truncation errors of order l/m^. On the downside, 
lattice spacings ai,a2 have to be considerably smaller and/or the masses mh in the numerical 
extrapolation mh — > mb have to be significantly smaller than the physical one. However, the 
combination Ljmh is always large compared to one and therefore a smooth extrapolation of the 
functions ai in l/(Ljmh) is expected and found numerically. A most promising approach is to 
combine the advantages of the two methods'^. In fact, let us compare 

"^b^("^b) = 4.347(48)GeV HQEtISSI 

quenched! 

m^^ifrib) = 4.421(67)GeV combinationSOl 

The comparison is meaningful because the experimental inputs (which matter in the quenched 
approximation) have been chosen the same. Let us finally note that the dominating uncertainty 
here is the one in the renormalization factor for the quark mass in the relativistic theory^IESl 
this can be improved. 

We emphasize that dynamical fermions are absolutely needed for a comparison to phe- 
nomenology. Presently, these strategies are being applied to A'^f = 2 QCD (just up and down 
sea). An intermediate resultSS shown in Fig. [3] demonstrates the precision of the most difficult 
continuum extrapolations. 

4-4 The CLS strategy 

While the numerical effort for the physically small lattices (see Fig. [2]) is quite modest for 
today's computer s, the theoretical setup is to be chosen with care. In particular the Schrodinger 
functional is used^^^ESl where Dirichlet boundary conditions in time play a prominent role. This 

again calls for a simple, very local action. 

Such an action, the 0(a)-improved Wilson actionSSEZHSl jg i^j^en also needed in the large 
volume part. The strategy of the CLS effort^^ is to use it and through its simplicity be able 
to profit from developments in algorithms'!^. This helps to diminish systematic and statistical 
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Figure 3: Example of a contin uum limit limag^o $*^*^''(ao, ii, mb) for a whole range of mb (left) and the resulting 
mass dependence on the right The variable z is the dimensionless RGI quark mass z = Mt,L. Curves are a 
low order HQET-representation. For details see the reference. 



errors by working at small lattice spacings and on physically large volumes. Furthermore per- 
turbative errors and in particular divergences (such as eq. ^) are removed by non-perturbative 
calculations in the Schrodinger functional . 

The present goal are simulations with A'f = 2 dynamical quarks, lattice spacings a = 
0.08, 0.06, 0.04 fm, sizes L = 2fm to 4fm and pion masses down to m.^ = 200 MeV. Con- 
cerning heavy flavours both the B-physics programme emphasized above and charm physics^^ 
are being investigated. We note that a considerable problem remaining will have to be solved 
along the way: the slow evolution of topological modes ^''^ in standard MC algorithms. 



5 Conclusions 

It is now widely recognized that lattice QCD is an important tool for computing phenomeno- 
logically relevant quantities. I have often heard complaints that progress is slow, but one has 
to recognize that the rich dynamics and the many relevant scales in QCD represent quite a 
challenge (Sect. [3]). Despite this we appear to approach the era of controlled precision predic- 
tions in time for an accurate analysis of present and future flavour physics data and models. In 
particular present worries about the use of the rooting trick and uncancelled power divergences 
can become history. 
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